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Abstract 


This  paper  considers  the  Lehmer  matrix  and  its  recursive 
analogue.  The  determinant  of  Lehmer  matrix  is  derived  ex¬ 
plicitly  by  both  its  LU  and  Cholesky  factorizations.  We  fur¬ 
ther  define  a  generalized  Lehmer  matrix  with  (i,j)  entries 
min  {iti+i,  Uj+i} 


9ij  = 


where  un  is  the  ?ith  term  of  a  binary 


max  {ui+i,Uj+i} 
sequence  {un} .  We  derive  both  the  LU  and  Cholesky  factor¬ 
izations  of  this  analogous  matrix  and  we  precisely  compute  the 
determinant. 


1  Introduction 

D.H.  Lehmer  (see  [2])  constructed  an  n  x  n  symmetric  matrix  A  = 
( aij)ij  whose  (i,j)  entry  is 

a  =  min  {i ,  j}  =  f  ifj  j  >  i, 

%3  max  {i,j}  \  j/i  i>  j- 

Define  the  second  order  recurrence  {Un  (p.  q)}  as  follows: 

Un  (P,  q)  =  pUn- 1  (p,  q )  -  qUn-2  {p,  q) , 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

2010 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-2010  to  00-00-2010 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


4.  TITLE  AND  SUBTITLE 

The  Lehmer  matrix  and  its  recursive  analogue 


6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  8.  PERFORMING  ORGANIZATION 

Naval  Postgraduate  School, Department  of  Applied  report  number 

Mathematics, Monterey, CA, 93943 

9.  SPONSORING/MONITORING  AGENCY  NAME(S )  AND  ADDRESS(ES )  10.  SPONSOR/MONITOR' S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

J.  of  Combinat.  Math,  and  Combinat.Computing  74  (2010),  193-205. 

14.  ABSTRACT 

This  paper  considers  the  Lehmer  matrix  and  its  recursive  analogue.  The  determinant  of  Lehmer  matrix  is 
derived  ex-  plicitly  by  both  its  LU  and  Cholesky  factorizations.  We  fur-  ther  de  ne  a  generalized  Lehmer 
matrix  with  (i;  j)  entries  gij  =  min  fui+1;  uj+lg  max  fui+1;  uj+lg  where  un  is  the  nth  term  of  a  binary 
sequence  fung  :  We  derive  both  the  LU  and  Cholesky  factor-  izations  of  this  analogous  matrix  and  we 
precisely  compute  the  determinant. 

15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

13 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


where  Uq  ( p ,  q)  =  0  and  U \  (p.  q)  =  1  for  n  >  1. 

As  an  interesting  example,  we  mention  that  the  set  of  natural 
numbers  can  be  obtained  from  the  sequence  {Un(p,q)}  by  taking 
p  =  2,  q  =  1.  Throughout  this  paper,  we  consider  the  case  q  =  —  1 
and  we  denote  un  =  Un  (p,  —1). 

We  now  define  an  n  x  n  generalized  Lehrner  matrix,  namely  Tn  = 
(9ij)i<ij<n  defined  below: 


mm{ui+i,Uj+i} 

max{ui+i,uj+i} 


where  un  is  the  nth  term  of  the  sequence  {un}  .  In  this  paper,  we 
obtain  the  general  LU  factorization  and  other  explicit  formulas  for 
both  the  Lehrner  matrix  and  its  recursive  analogue. 

The  Lehrner  matrix  is  part  of  a  family  of  matrices  known  as 
test  matrices,  which  are  used  to  evaluate  the  accuracy  of  matrix 
inversion  programs  since  the  exact  inverses  are  known  (see  [1,2]).  It 
is  hoped  that  our  generalized  Lehrner  matrix  will  add  to  the  literature 
of  special  matrices  with  known  inverse. 


Uj+l 
Uj  + 1 

ui+ 1 


if  3  >  i, 


if  i  >  j. 


2  The  Lehrner  Matrix 

We  start  by  obtaining  the  LU  factorization  of  the  Lehrner  matrix 
A.  Using  the  inverses  of  L  and  U,  we  obtain  the  explicit  form  for 
the  inverse  of  A,  whose  inverse  is  well-known,  thus  obtaining  another 
proof  of  this  result. 

We  define  the  n  x  n  invertible  lower  triangular  matrix  L  =  (£tj) 
where  £tj  =  j/i  for  i  >  j  and  0  otherwise.  Next,  we  define  the  n  x  n 
invertible  upper  triangular  matrix  U  =  (uij)  with  uVJ  =  for 


2 


i  <  j  and  0  otherwise.  For  example,  when  n  =  5,  we  get 
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The  following  result  holds. 

Theorem  1.  For  n  >  0,  the  LU  factorization  of  Lehmer  matrix  is 
given  by 

A  =  LU 

where  L  and  U  were  defined  previously. 


Proof.  We  split  the  proof  into  three  cases. 
Case  1 :  i  =  j.  By  Y?k=i  (2fc  —  1)  =  t2,  then 


O'ii  — 


-  Et2^  -  E  ^ 


fc=l  Al=l 

Case  2:  i  >  j.  Thus 

n  j 


,'2 


k= 1 


k  (2k  —  1) 


aij  —  ^2^ikUkj  —  ^2^ikUkj  ~  ■  j,  ■ 

k= 1  /c=l  fc=l 


-  E 

k~- 

Case  3:  j  >  i.  Then 


3  2k  -  1  1  „  j 


%■)  ii 

k= 1  J  J  fc=l 


-V2fc-1 
)  1  ^ 


ai]  ~ 


^  '  ^ik  ^kj  ^  '  I'ik  ^kj  '  . 


fc=l 


fe=l 


fc=l 


k  (2k  —  1) 
i  kj 


V  — u-i 

Z — ^  i  I  ]  i  Z — J 


11  11 

k= 1  J  fc=l 
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which  completes  the  proof.  □ 

We  display  an  example  below: 
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As  a  consequence  of  Theorem  1,  we  obtain  an  explicit  value  of 
the  determinant  of  the  Lehrner  matrix  in  the  following  corollary. 


Corollary  1.  For  n  >  0, 


det  A 


(2  n)\ 
2n  (n!)3 


Proof.  The  proof  follows  from  the  LU  factorization  of  matrix  A  by 
considering  det  A  =  det  U  =  n  7=  i  □ 


by 


The  nth  Catalan  number  is  given  in  terms  of  binomial  coefficients 

c  1  (1  n\  =  (2  n)! 

n  n  +  l\nj  (n+l)!n!" 


Thus  we  may  note  that 


det  A 


(n  +  1) 
2nn! 


Cr, 


We  continue  our  analysis  by  determining  the  L\Lj  (named  after 
Cholesky)  factorization  of  the  Lehrner  matrix,  where  L\  is  a  lower 
triangular  matrix.  The  Cholesky  factorization  was  obtained  for  a 
different  kind  of  matrix  defined  using  binary  sequences  by  the  second 
author  in  [3]. 


4 


Theorem  2.  The  Cholesky  factorization  of  the  Lehmer  matrix  is 
given  by 


A  =  L  iLf 


where  L\  =  (f%j)  is  a  lower  triangular  matrix  with  = 
all  i  >  j. 

Proof.  If  i  >  j,  then 


_  VWfA 


for 


t  Mr  =  t  Mr  =  t  VW^~1 

r= 1  r= 1  r=  1 


j 


rP{2r-1] 

J  r= 1 


If  z  =  j,  then 


Y.A=Y.fl  =  E 


r=l 


r=l 


-^(2r-i)  =  ^  =  I- 


r=l 


r=l 
2 


Finally,  if  i  <  j ,  then 


Clij  — 


^  ]  firfjr  —  firfjr  —  .  .  (2r  1) 


r=l  r=l 

which  proves  the  theorem. 


r=l 


□ 


As  an  example,  for  n  =  5  and  p  =  1  (the  Fibonacci  sequence 
case),  we  have 
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By  Theorem  2,  we  find  that,  since 
1 2  _  rrn  2i— 1  _  (2w)! 


det  (A)  =  nr=i  /»  =  Ylt=i 


2ri(n!) 


7; 


A  =  L\Lj ,  we  have  that 
that  is,  Corollary  1. 


3  The  Inverse  of  the  Lehmer  Matrix 


Now  we  find  an  explicit  formula  for  the  inverse  of  the  Lehmer  matrix. 
For  this  purpose,  we  use  its  LU  factorization  as  A -1  =  U~1L~1.  We 
first  derive  the  inverses  of  the  matrices  L  and  U. 


Lemma  1.  Let  L  1  =  (ttj)  denote  the  inverse  of  L.  Then 


f  1  if  i  =  j, 
Uj=l  -f  ifi  =  j  + 1, 

[  0  otherwise, 


Proof.  The  proof  can  be  easily  checked  from  the  product  L  1L.  □ 

Lemma  2.  Let  U~l  =  (w%])  denote  the  inverse  of  U.  Then 


2i—l 
i(i+ 1) 
2i+l 
0 


if  i  =  3 
ifi  +  l  =  j, 
otherwise, 


Proof.  The  proof  follows  from  the  product  U  lU. 


□ 


The  inverse  of  the  Lehmer  matrix  is  found  in  the  following  theo¬ 
rem. 


Theorem  3.  For  n  >  0,  let  A  1  =  (bij),  then 


bij 


4  i3 

4i2  — 1 
n 2 

2n-l 

i(*+l) 

2i+l 

0 


if  i  =  j  <  n 
ifi  =  j  =  n, 

if  |  i  -  j  |  =  1, 
otherwise, 
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Proof.  Since  A  1  =  U  lL  1,  using  the  previous  two  lemmas,  we 
obtain  for  1  <  i  <  n  —  1, 

n 

ba  =  ^  '  ^iktki  =  'Wii  “I- 

fc=l 

i2  *  (i  + 1)  i  i2  i2  4i3 

2i  -  1  +  2i  +  l  (i+1)  ~~  2 i  -  1  +  2i+  1  ~~  4i2  -  1 ' 

2 

When  i  =  j  =  n,  it  is  easy  to  see  that  bnn  =  wnn  =  2n-i  •  ^ 
i  =  j  +  1 ,  then 

n 

bi-\-\,i  —  ^  '  'UJi+l.ktki  — 

fc=l 

(i  +  1)2  /-i\_  i(i  +  l) 

2i  +  1  \i  +  1  y  2i  +  1 

The  last  case  j  =  i  +  1  can  be  similarly  done,  and  the  proof  is 
complete.  □ 

Therefore  we  recover  the  known  fact  that  the  inverse  of  the  Lehrner 
matrix  is  a  symmetric  tridiagonal  matrix. 

We  give  the  following  example  as  a  consequence  of  the  above 
theorem:  for  n  =  4, 


7 


We  also  give  a  relation  between  the  terms  of  inverse  of  the  Lehrner 
matrix  and  triangular  numbers.  Recall  that  the  nth  triangular  num¬ 
ber  Tn  is  defined  as  the  sum  of  the  first  n  natural  numbers,  that 
is,  Tn  =  1  +  2  +  •  •  •  +  n  =  id+ti).  We  can  re-write  A~l  =  (bjj)  as 

bij  =  -  Jqfr  for  | i  -  j\  =  1,  and  bu  =  . 


4  Recursive  Analogue  of  the  Lehmer  Matrix 


In  this  section  we  investigate  the  same  questions  for  our  general¬ 
ized  recursive  analogue  of  the  Lehmer  matrix  Tn  defined  in  the  first 
section,  namely,  J~n  =  {[hj)'- 


min  {uj+i,Uj+i} 
max{«i+i,uJ+i} 


Uj+ 1 

Uj  + 1 
Uj  + 1 

Ui- 1-1 


if  j  >  h 


if  i  >  j. 


where  un  is  the  nth  term  of  the  sequence  {un} 
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following  form: 
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In  order  to  give  the  LU  factorization  of  the  matrix  Tn ,  we  define 
two  triangular  matrices. 

Define  the  n  x  n  unit  lower  triangular  matrix  L2  =  (cjj)  with 
ca  =  for  all  i  >  j  and  =  0  for  all  i  <  j. 

J  Ui-\- 1  J 


For  example,  when  n  =  5,  the  matrix  takes  the  form: 
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Before  defining  an  upper  triangular  matrix  for  the  LU  factorization 
of  the  matrix  J~n.  we  need  to  introduce  a  new  sequence  {tn}  by  the 
following  relation: 


tn  =  (p-  1)  Un  +  un- 1,  that  is,  tn  =  un+ 1  -  un,  n  >  1, 


where  u„  is  defined  as  before. 


u  2 

Uj  + 1 


Define  the  n  x  n  upper  triangular  matrix  U2  =  ( dij )  with  d\j  = 
for  1  <  j  <  n,  dij  =  for  1  <  i  <  j  <  n. 


■i+luj+l 


From  the  definition  of  the  sequence  {tn}  ,  we  rewrite  the  matrix 

for  1  <  i  <  j  <  n. 


U2  with  d\j  =  for  1  <  j  <  n,  dij  = 


UJ  1  1  —uj 


Uj  + 1  ”  —  J  —  *•?  Wj+lWj  +  1 

For  example,  when  n  =  4,  the  matrix  takes  the  form: 
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Theorem  4.  For  n  >  0,  the  factorization  of  matrix  Tn 
given  by 

tFn  =  L2U2, 


{dij)  ^ 
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where  U2  and  L2  were  defined  previously. 


Proof.  Let  L2U2  =  ( hij )  .We  consider  two  cases,  i  >  j  and  i  <  j.  For 
the  first  case,  we  write 


hij  —  E  C-im dmj  —  ^  ^  Cimdmj 


m=  1 


m=l 


c  d]_-  +  f Um+i  u™) 

lj  m=2  \  Ui+1  Um+lUi+l 


Uo 
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1 


^i+l^j+1  V'i+l'U'j+l  _ 


^  ^  (fim+ 1  ^m) 


m= 2 
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/  2  2\  uj+l 

K+i  -U2)  =  — —  =  gij. 

Ui+1Uj+l  Ui+lUj+l  J T  tti+1 


If  i  E  j ,  then  similarly 

n  i 

hij  —  y  ^  Cimdmj  —  y  '  Cirnd 


mj 


m= 1 


m=l 


—  Cil 


+  X] 


/  Mm+l  (^m+l  Um) 


_2  y  ^i+l 


+ 


1 


uj+ 1 

and  the  claim  is  shown. 


'Uj+l'^j+l  Ui+lUj+l  m=2 

a  1  ■  1  _ 

—  ) 


X  (“m+1  “  “m) 


□ 


Now  we  can  find  the  value  of  det  (Fn)  by  considering  its  LC/ 
factorization. 


Corollary  2.  For  n  >  0, 

n 

det  (Tn)  =  J] 


i= 2 


u: 


i+ 1 
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As  a  special  cases  of  the  matrix  Tn ,  we  take  the  matrix  F°  ob¬ 
tained  using  the  Fibonacci  sequence,  that  is,  Fn+\  =  Fn  +  Fn-i,  Fq  = 
0,  F\  =  1.  The  determinant  of  this  matrix  becomes 


det  (J%) 


Fn-i-Fn+2! 

2(Fn+1!f  ’ 


where  Fn\  is  the  Fibonomial  factorial,  that  is,  Fn\  =  F1F2  ■  ■  ■  Fn. 

Next  we  give  the  Cholesky  factorization  of  the  generalized  Lehmer 
matrix  Fn.  For  this  purpose  we  define  a  lower  triangular  matrix  L3  = 

(mu)  with  mu  =  for  1  <  i  <  n,  mu  =  —^—\U2- —  u2  for 
V  iJJ  ui+l  ——11]  u.+1  y  j+1  j 

1  <  j  <  i  <  n  and  0  otherwise. 

When  n  =  4,  the  matrix  L3  takes  the  form: 
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The  proof  of  the  next  theorem  is  analogous  to  the  proof  of  The¬ 
orem  4,  so  it  will  be  omitted. 

Theorem  5.  The  Cholesky  factorization  of  the  recursive  analogue 
of  the  Lehmer  matrix  is  given  by 

Fn  =  A3L3 

where  L3  is  the  lower  triangular  matrix  defined  previously. 


5  The  Inverse  of  the  Generalized  Lehmer  Ma¬ 
trix 

Here  we  give  the  inverse  of  the  recursive  analogue  of  the  Lehmer 
matrix  F~l  by  considering  its  LU  factorization.  Before  this,  we  give 
the  inverses  of  the  matrices  L2  and  U2  in  the  following  lemmas,  stated 
without  proofs,  as  they  are  immediate. 
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Lemma  3.  Let  U2  =  (r%)  denote  the  inverse  of  U-2 ■  Then 


10  ij  =  < 


if  i  =  j  =  l 


,j2i+ 1 


if  1  <  *  =  3, 

ui  ui+ 1 

«i+l“i+2  -r  ■  |  i  —  • 

“i+l  ai+2 


otherwise , 


Lemma  4.  Let  L2  =  (t^/)  denote  the  inverse  of  L.  Then 

i  */  * = j, 

4  =  <(  ifi  =  j  +  1> 

0  otherwise, 


Thus  the  inverse  of  the  matrix  Tn  is  found  in  the  following  the¬ 
orem. 

<4_ 

Uo—uZ 


Theorem  6.  For  n  >  0,  let  Tn  1  =  (qij),  then  qn  =  2_3  2  ?  qnn  — 


Un+1  /r.  .  „  _  ~  „  .  _  ^i+l^+2 

,.2  _,#2  5  W+l  —  </2+l,2  —  ,.2  _ 7/2 

^n+l  an  ai+l  ai+2 


for  1  <  i  <  n  -  1,  qu  = 


2  /  2  2\ 

M'+2w+2  ^  2 — y  for  2  <  i  <  n  —  1  and  0  otherwise. 


{ui+l-ui){ui+2-ui 


+  1) 


Proof.  Since  3~n  1  =  U2  1L2  1,  the  proof  follows  from  the  previous  two 


lemmas  and  from  matrix  multiplication. 


□ 


For  example,  for  n  =  4, 


f- 1  — 
""  5  — 


3  “2 
yyU3 


U2U3 


3  2  2 
uZ—uZ 


U3U4 
„.2  ,.2 


^3^4 


uif—u 


4  2  2 
ui—uZ 


U4U5 


0 

0 

W4U5 

““ 2 - 2 

uZ 


6  Further  comment 


With  a  bit  more  care,  one  can  certainly  remove  the  constraint  q  =  —  1 
on  the  sequence  Un,  and  prove  similar  results  like  in  the  present  paper 
for  the  corresponding  generalized  Lehnrer  matrix. 
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